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ABSTRACT 

We have observed pulsar wind nebulae in the three supernova remnants Gil. 2— 0.3, G16. 7+0.1, and 
G29.7— 0.3 at 89 GHz with the Berkeley-Illinois-Maryland Association Array, measuring total flux 
densities of two of them for comparison with archival data at other frequencies. In G16.7+0.1, we find 
a break in the spectrum of the PWN at ~ 26 GHz. In G29.7— 0.3, our data suggest a break in the 
integrated spectrum of the central nebula at ~ 55 GHz, lower than previously estimated. However, 
we have found spatial structure in the spectrum of this nebula. The emission to the north of pulsar 
J1846— 0258 has a broken spectrum, with break frequency < 100 GHz, consistent with a conventional 
pulsar-powered nebula. The emission to the south of the pulsar has a near-power-law spectrum from 
radio to X-rays: this component may be unrelated to the PWN, or may be evidence of asymmetries 
and/or time evolution in the pulsar's energy output. We present 89 GHz images of each remnant. 

Subject headings: pulsars: individual (J1846 — 0258) — radio continuum: ISM — SNR: individual 
(Gil. 2-0.3, G16. 7+0.1, G29. 7-0.3)— supernova remnants 



1. INTRODUCTION 

Pulsars are born spinning with up to 10 51 erg of ro- 
tational kinetic energy — as much energy as in the su- 
pernovae in which they are born. This vast reservoir of 
energy is ultimately deposited into the ambient medium 
through the pulsar's relativistic wind. In many cases, 
this interaction between the pulsar and its environment is 
directly observable, in the form of a synchrotron-emitting 
pulsar wind nebula (PWN) , of which the Crab Nebula is 
the best-known example. While many pulsars have their 
beams directed away from us, PWNe radiate isotropi- 
cally. Thus PWNe are powerful probes of pulsars and 
their energy loss, even when the pulsar itself cannot be 
detected. 

At radio frequencies, PWNe typically have flat spec- 
tra —0.3 < a < (S v oc v a ), but in the X-ray 
band we generally see a < — 1. It is therefore com- 
monly presumed that PWNe have at least one spec- 
tral break at intermediate wavelengths. Such breaks 
are expected from theoretical considerations, and re- 
sult from a combination of synchrotron losses and the 
t ime -evolution of the pulsar's changing energy out- 
put JPacini fc Salvatilfl973t IRevnolds fc Chevalierlll984t 
IWoltier et al.H1997lK Locating these breaks gives insight 
into the physical conditions of the pulsar wind: in PWNe 
powered by young (< 5 kyr) pulsars, the spectral break 
can be used to infer directly the nebular magnetic field 
stren gth (e.g. iManchester. Stavelev-Smith. fc Kestevenl 
1993), while spatial variations in the break energy can be 
used to identify and map out the processes of particle dif- 
fusion and radiative losses within th e flow iJAmato et alJ 
2000; [Bock. Wright, fc Dickell200l . 
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Most flux density measurements of PWNe to date have 
been in the X-ray and low-frequency radio bands. To pin 
down the break frequencies it is necessary to measure flux 
densities in the part of the spectrum near the spectral 
breaks, i.e. at millimeter and infrared wavelengths. 

Here we report on observations of three PWNe for 
which extant data at radio and X-ray wavelengths 
implied a spectral break in or below the mm band. 
The PWNe, in the supernova remnants Gil. 2— 0.3, 
G16. 7+0.1, and G29.7-0. 3 (Kes t even 75) , were iden- 
tified from the catalog of iGreenl (|2004fl as being suffi- 
ciently bright and compact to allow imaging with mil- 
limeter interferometers. Two of the PWNe (in Gil. 2— 0.3 
and G29.7— 0.3) contain you ng X-ray pulsars (jTorii et alJ 
IT9971 IGot^helf^t^l^. Only a few PWNe have 
been reliably imaged at millimeter wavelengths with sin- 
gle dish tel escopes, and prio r to this study only one 
fG21. 5-0.9: IBock et aTll2001J) had been imaged with a 
millimeter interferometer. 

We have been able reliably to image two of the sources 
observed (in the supernova remnants G16. 7+0.1 and 
G29.7— 0.3), obtaining their total flux densities at 89 
GHz. This has allowed us to determine the frequency of 
their spectral breaks at radio wavelengths using archival 
data and assuming the presence of a solitary break. We 
have been able to resolve the PWN in G29.7— 0.3, letting 
us compare the millimeter-wave structure of the PWN 
with that in the X-ray band and at lower radio frequen- 
cies. We also present a millimeter- wave image of the SNR 
Gil. 2— 0.3, which has a PWN with radio spectr al index 
a = -0.25lg;5o llTam. Roberts, fc Kasoil l200l . How- 
ever, its extension compared to the spatial sensitivity of 
the array precludes any quantitative analysis. 
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2. OBSERVATIONS 

Observations were obtained a t 88.6 GHz using the 
BIMA Array 2 (|Welch et al.lll99l| during 2002 and 2003 
(Table QJ. To achieve suitable uv coverage, we used 
the technique of multi-frequency synthesis over both 800 
MHz sidebands of the local oscillator. We m easured only 
the left circular polarization (in the sense of lIEEFj|l969j) 
in the assumption that circularly polarized emission from 
the sources would be negligible. 3 In order to image fully 
the extended structure, and to assist with the recovery 
of low spatial frequency information, we imaged each 
source with a 7-point hexagonal mosaic. The primary 
beam of the 6.1-m antennas is Gaussian with FWHM 
2f 13. The short interferometer spacings (as low as 6.1 m 
when antennas are nearly shadowed) available from the 
array provide significant sensitivity on spatial scales up 
to an arcminute. The shortest physical baselines are 8.1 
m. 

The bright quasars QSO B1730-130 (for Gil. 2-0.3 
and G16. 7+0.1) and QSO B1741-038 (for G29. 7-0.3) 
were observed as phase calibrators at 25 to 30 minute 
intervals, after 3 complete cycles of 7 pointings. We ex- 
pect the astrometric accuracy of the images to be bet- 
ter than 1 arcsec, a precision routinely obtai ned with 
the BIMA Array l)Loonev fc Hardcastlel 120001) . Single 
sideband system temperatures during the observations 
were between 200 and 600 K (scaled to outside the 
atmosphere). We used the on-line absolute flux den- 
sity scale determined from many observations of plan- 
ets, which have been found on average to be stable over 
time l|Muhleman fc Bergell99l|) . As a check we observed 
Uranus for 10 minutes during each observing session, 
finding a day-to-day scatter for the mean gain on all an- 
tennas between 0.95 and 1.12, implying multiplicative 
corrections to the flux density scale of between 0.91 and 
1.25. The mean correction would be 1.09 for observa- 
tions of G16.7+0.1 and 1.05 for G29. 7-0.3. The gains 
for each sideband were consistent within a few percent 
on each day. In this paper we retain the on-line scale, 
and note that with several days' observations contribut- 
ing to each image our measurements are consistent with 
an uncertainty in the overall flux density scale of about 
10%. 

Im aging was done with the MIR IAD software pack- 
age (iSault. T euben. fc Wright! 119951). u sing SDI decon- 
volution l|Steer. Dewdnev. fc Itohl 1 19841k such CLEAN- 
based algorithms are generally superior to maximum en- 
tropy methods for low signal to noise ratios ijHelfer et alJ 
2003) . The multiple pointing centers were combined into 
a mosaiced image corrected for the primary beam. To 
assess the reliability of flux density recovery when imag- 
ing these sources, simulations of each observation were 
made with MIRIAD, using the actual observations as a 
template. The results of the simulations are presented 
with the data on each source in the next section. 

3. RESULTS 

2 The BIMA Array was operated until June 2004 by the 
Berkeley-Illinois-Maryland Association with support from the Na- 
tional Science Foundation. 

3 The BIMA Array natively measured one linear polarization. 
To convert to circular polarization, 90-GHz quarter-wave plates 
were used. 




18 h ll m 40 s 36 s 32 s 28 s 24 s 20 s 16 s 
RA (J2000) 

Fig. 1.— Image of Gll.2-0 j at 89 GHz (contours) and 8.6 GHz 
(grayscale; IRoberts et al. 2003). The contours are at 10, 20, and 
30 mjy beam" 1 . The synthesized beam at 89 GHz is shown at the 
lower right. 



3.1. Gil. 2-0.3 

The BIMA image of Gil. 2— 0.3 is shown with con- 
tours in Figure [H o verlaid upon an 8.6 GHz image of 
IRoberts et al.l l)2003jh Based on positional coincidence 
with the 8.6 GHz emission, we appear to have detected 
the PWN in this remnant (the single contour at the cen- 
ter of the image). However, the extended structure in 
the PWN is not well imaged by the array: the model- 
ing described in section[2]showed that we have measured 
only approximately 10% of the emission, indicating that 
the PWN is nearly resolved out on all baselines. Given 
the magnitude and unreliability of the correction factor 
required to estimate the true flux density, we do not be- 
lieve that these data can provide a reliable estimate of 
the flux density of the PWN. High sensitivity single dish 
measurements would be required to investigate this ob- 
ject further. 

3.2. G16. 7+0.1 

The PWN in this remnant (Figure |5J) is detected as 
a peak of 14 mjy beam" 1 at 8 h 20 m 57 s 2 -14°19'34" 
(J2000). The feature is extended: an elliptical Gaussian 
fitted to the feature has a size of 51" x 20" (at position 
angle 6° east of north), and an integrated flux density of 
27 mjy. The siz e and position are c onsistent with the 
measurements of iHelfand et alJ l)1989[L 

To assess the reliability of this measurement we have 
modeled the response of the array to two simple sources: 
a 25" x 7" Gaussian with major axis 10° west of north 
and a 45" x 15" Gaussian with major axis 20° west of 
north, corresponding approxim ately to the 50% an d 20% 
contours in the 2 cm image of IHelfand et alJ 1)19891 ). In 
deconvolved images we recovered 97% and 80% of the 
total emission respectively, without taking into account 
a small residual negative background. We conclude that 
we should have measured at least 90% of the emission 
from this object, and consider that we have reliably de- 
termined the flux density of the PWN at this frequency. 
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TABLE 1 

Summary of observations 



SNR 


Antenna 
Config. 


Date(s) 


Time on 
Source 
(hr) 


Image rms 
(mjy beam -1 ) 


Synthesized Beam 
FWHM size PA 

(") (°) 


Gil. 2-0.3 


D 


2002 Jun 10 


3.7 


< 10 a 


18.9x17.6 19.0 


G16.7+0.1 


D 


2002 Jun 2, Sep 10; 2003 Sep 1,4 


14.3 


~ 4 


25.2x18.7 9.7 


G29.7-0.3 


D 


2002 May 30, 2003 Aug 31 


8.2 


1.7 b 


13.3xl0.4 b 17.1 b 




C 


2002 May 9,16 


7.9 








B 


2002 Mar 15 


5.6 







a Low-level emission throughout the field makes the rms difficult to determine accurately 
b For C and D configuration observations 
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Fig. 2.— Image of G16. 7+0.1 at 89 GHz. The contours are at 4, 
8, and 12 mjy beam -1 . The synthesized beam at 89 GHz is shown 
at the lower right. 



Taking into account the uncertainty in the flux density 
scale (@ and the low signal to noise we estimate the to- 
tal uncertainty in our measurement to be no more than 
20%. 

3.3. G29. 7-0.3 (Kesteven 75) 

The BIMA image of G29.7-0.3 is presented in 
Figure EH For comparison, the 1.4 GHz image of 
IHelfand Collins, fc GottheU pilK ihi is shown also. This 
object is sufficiently bright at millimeter wavelengths 
that we have been able to image it using both the BIMA 
C and D configurations. However, the data obtained 
with the B configuration did not add significantly to the 
results: they have not been included in the image in Fig- 
ure |3 nor in the quantitative analysis that follows. In 
the figure, both the PWN (at the image center) and the 
partial shell (to the south and east) can easily be seen 
at each frequency. As expected, the PWN is relatively 
brighter than the shell at the higher frequency, i.e. it has 
a flatter spectral index. Note th at the PWN is centra l 
in the corresponding X-ray image ijHelfand et al] l2003b). 
which also shows the northern extent of the SNR's shell. 

We have measured an 89 GHz flux density for the PWN 
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Fig. 3.— Images of G29.7-0.3 at 1.4 GHz IHelfand et alJ2003ri 

contours 2, 5, 10, 20, 50, 80 mjy beam -1 on logarithmic grayscale) 
and 89 GHz (linear scale; contours at 2, 5, 10, 20 mjy beam - 1 ). 
The position of the pulsar J1846-0258 IGotthelf etalllSoOCh is 
marked with a cross. The synthesized beams are shown at the 
lower right of each image. 
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of 80 mjy in 0.6 arcmin 2 (within approximately the "0 
mJy beam -1 contour" ), taking into account a local back- 
ground (a "negative bowl" ) within the SNR shell of about 
—36 mJy arcmin -2 . To quantify the degree to which the 
array is sensitive to the extended emission, we simulated 
an observation of the PWN using the 1.4 GHz image 
(IHelfand et al.l l2003b) as a model. 4 We added Gaussian 
noise to the undeconvolved image to give a similar signal 
to noise ratio. After deconvolution, we recovered 92.4% 
of the flux density in the simulated observation. Ap- 
plying the correction for interferometric filtering of the 
extended emission thus determined, we obtain a flux den- 
sity for the PWN of 86 mJy. In the simulation the neg- 
ative background was much less pronounced than in the 
BIMA observation: this artifact is thus probably due pri- 
marily to poor deconvolution of emission from the SNR 
shell, which was not included in our simulation. We con- 
sider an upper limit to the imaging errors to be 10%; 
combining this with the 10% flux calibration uncertainty 
leads to a total uncertainty of 15%. 
At the nearby frequ ency of 84 GHz, the observations 
of iSalter et alJ l)1989ft imply an integrated flux density 
of 159 mJy, inconsistent with the present data. How- 
ever, the earlier observations were made using a single 
dish with a beam of 76" x 70", substantiall y larger than 
the characteristic source size (~ 30" x 20"). ISalter et all 
corrected their measurement for the smaller source size 
to derive their integrated flux density. However, it seems 
possible that Salter et al. also measured some shell emis- 
sion. The advantage of the present interferometric mea- 
surement is that it would have filtered out any contam- 
inating emission from the SNR's shell, and should pro- 
vide a more reliable determination of the flux density 
attributable to the PWN alone. 

4. DISCUSSION 
4.1. G16.7+0.1 

A spectrum of the PWN in G16. 7+0.1 is shown in Fig- 
ure 0] There are few measurements available of the flux 
density of this PWN. The solid line at low frequencies 
represents the radio sp ectrum derived from observations 
at 6, 20, and 90 c m bv IHelfand et alJ (119891). Recent X- 
ray measurements (Hclf and. Agiieros. fc Gotth clf 2003a) 
imply that there must be at least one break between the 
radio and X-ray regimes; the present measurement im- 
plies that this is at about 26 GHz (central dotted line). 
In the figure, the gray shaded area represents the region 
of intersection of the Helfand et al. radio spectrum and 
a simple radio/X-ray spectrum (i.e. assuming a single 
break between X-ray and radio wavelength) consistent 
with the 89 GHz and X-ray uncertainties. This pro- 
vides one way to quantify the uncertainty in the break 
frequency: the region has minimum and maximum fre- 
quencies (shown by the outer dotted lines) of 16 and 42 
GHz respectively, leading to an estimate of the break 
frequency of 26±^ GHz. 

The PWN in G 16. 7+0.1 has one of the lowest known 
break frequencies in a PWN: breaks at lower or compa- 
rable frequencies have been reported only in G27.8+0.6 
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Fig. 4. — Integrated spectra of the pulsar wind nebulae in 
G16. 7+0.1 and G29.7— 0.3. In regions where the spectra are well 
determined they are shown with unbroken lines. Elsewhere simple 
straight line spectra implied by sparse measurements are shown 
by a dashed line. Uncertainties are shown where available; at X- 
ray wavelengths these are indicated by black shading. The regions 
within which spectral breaks are implied are shown by gray shad- 
ing, and bounded by the vertical dotted lines (see text for details). 



(between 5 and 10 GHz; iReich. Furst. & Sofud Il98l 
and G 74.9+1.2 (between 10 and 30 GHz: lMorsi k. Reicbl 
l!987a|L Surprisingly, the spectrum within the X-ray 
band (a = —0.17 ± 0.29) is flatter than the radio to 
X-ray spectrum. Presumably this is not the result of 
inverse Compton scattering, which is generally expected 
to be significant at 7-ray wavelengths. This phenomenon 
has been seen in one other PWN, t hat around the pulsar 
PSR B1757-24 (jKaspi et al.l2001() . 5 In that case, the X- 
ray and radio spectra are also approximately flat, with a 
steeper spectrum required to con nect the radio and X-ray 
observations. However, whereas iKaspi et al.l found that 
the X-ray spectrum was just marginally consistent with 
the implied radio to X-ray average spectrum, this possi- 
bility is very unlikely in the case of G16. 7+0.1, where the 
average radio to X-ray spectrum is not consistent with 



4 We chose this image for a model owing to its higher signal 
to noise. Recalling the fiat PWN spectrum in the radio, we do 
not believe that the difference between the PWN shape at the two 
frequencies will significantly affect our result. 



5 A step like that in the spectrum of the Crab Nebula at in- 
frared/o ptical wavelengths, thought d ue to dust reradiation or ex- 



tinction (Manchester & Tavlor 1977), would not show up in the 
bulk radio to X-ray spectral index in the same way. 
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the X-ray data even w ithin the 90% uncertainty limits of 
IHelfand et a,l]Ef&. 

One explanation for spectra of this shape is given 
bv iRevnolds fe Chevalieif |l984'l. They calculate that 
a combination of older particles suffering synchrotron 
losses due to compression of the magnetic field by the 
reverse shock and particles injected after the shock can 
give rise to a spectrum of this shape. However, whereas 
the re serve shock is expecte d to arrive after about 10 4 
years, Hclf and et all {2Q£)3a) have estimated G16.7+0.1 
to be only about 2000 years old. Therefore this seems 
like an unlikely explanation for this source. 

4.2. G29.7-0.3: Integrated Spectrum of the PWN 

An integrated spectrum of the central nebula in 
G29.7— 0.3 is plotted in FigureJH This object was first 
seen bv lBecker fc Kundul l)1976|) . who measured a radio 
spectral index of —0.27, flatter than that of the SNR 
shell. VLA observations ijBecker fc Helfandlll984(l con- 
firmed the discovery but implie d a flat (a = 0.0) spec- 
trum. Hunt et al. (as reported by Sal ter et al.ll989D from 
different VLA observations obtained integrated flux den- 
sity measurements that imply a radio spectral index of 
—0.30. Both sets of VLA observations are plotted in Fig- 
ure 0] We have measured a flux density of 343 mJy for 
the PW N in the recent 1.4 GHz data of IHelfand et all 
l)2003b(l . consistent with t he Salter et al. r esult. 

At higher frequencies, iMorsi fc Reicbl l)1987br) mea- 
sured 0.28 Jy from the PWN at 30 GHz, but note that 
they could have suffered from confusion (their observa- 
tions had resolution 26'.'5), leading to an overestimated 
result. It seems possible that they measured some of the 
broad underlying shell-related emission as well, which 
was not imaged at 89 GHz, and could have been poorly 
imaged at lower radio frequencies. This could have led 
to a furt her overestim ate. As noted earlier f i)3.3|l we con- 
sider the [Salter et al. 84 GHz measurement also to be an 
overestimate, for similar reasons. 

Combined with the Salter et al. lower- frequency mea- 
surements, the X-ray data of Hclfan d~et alJ l)2003b|) . and 
assuming a simple unbroken spectrum between radio and 
X-ray wavelengths, our observations imply a spectral 
break at approximately 55 GHz (shown by the right-hand 
vertical dotted line in the figure). Assuming the break 
to be sharp, the integrated flux density of the PWN at 
this frequency would be about 120 mJy beam^ 1 . Allow- 
ing the spectrum above and below the putative break to 
vary within the uncertainties of the measurements im- 
plies that the break should fall within the shaded gray 
area, with a minimum frequency of 21 GHz (left-hand 
vertical line). The break is not strictly constrained at 
high frequencies, since a straight line is just consistent 
with the low-frequency Hunt et al. measurements and 
our 89 GHz measurement. However, it is likely that the 
break is below 89 GHz. 

We c ould alternatively have us ed the centimeter-wave 
data of iBecker fc Helfandl l|1984D and obtained a some- 
what lower break frequency. However, gi ven the agree- 
ment b etween the data of Salter et al. and IHelfand et al.l 
(2003b), and the fact that the Salter et al. data were 
quoted with uncertainties, we prefer the later measure- 
ments. 

4.3. G29.7-0.3: Lobes of the PWN? 



The high resolution X-ray data clearly show two lobes 
of emission, to the north and south of the pulsar. The 
northern lobe is brighter at X-ray wavelengths. A com- 
parison of images of the PWN at 1.4 GHz, 89 GHz and 
in X-rays (Figure 0) reveals that the PWN spectrum is 
also asymmetric. At 1.4 GHz, the PWN is brighter south 
of the pulsar, while at 89 GHz and above it is brighter 
north of the pulsar. These two components may have 
different origins, causing an integrated spectrum to be 
difficult to interpret. 

In Table [2J we show flux densities of these northern 
and southern lobes, integrated over the regions shown 
in Figure The regions chosen were a compromise 
between the competing desires to minimize contamina- 
tion from the background and to include the majority 
of emission at each frequency. The area around the pul- 
sar was excluded from all datasets, avoiding the need to 
estimate an X-ray background for normalizing the inte- 
grated areas. A t 89 GHz the same ba ckground as in 
was subtracted. IHelfand et alJ l)2003bl) were not able to 
find a significant difference between the X-ray spectra of 
the northern and southern lobes, so the X-ray flux den- 
sities quoted correspond to the luminosity reported by 
IHelfand et all apportioned according to the total counts 
in the two regions, assuming a uniform X-ray spectrum 
across the PWN. The 89 GHz uncertainties quoted are as 
in sectional those for X-rays are the standard errors for 
a Poisson distribution, and we assign an uncertainty of 
10% to the 1.4 GHz measurements, to account for imag- 
ing errors. This is probably an overestimate of the true 
uncertainty at 1.4 GHz. 

Given the differing resolutions and systematic errors 
of the three datasets 6 , the results are necessarily quali- 
tative. However, as may be seen in Figure El the north- 
ern lobe not only has a flatter spectrum on average, but 
shows strong evidence for break. Meanwhile, the spec- 
trum of the southern lobe may be unbroken. In the fig- 
ure, the uncertainties at radio frequencies are of order 
the point size, while those of the X-ray values are much 
smaller. The excellent positional agreement of the SNR 
shell at 1.4 and 89 GHz confirms the registration of the 
images. We believe that the difference in spectral shape 
between the lobes is reliable. 

The same result may be obtained qualitatively by con- 
volving the 1.4 GHz and X-ray observations to the reso- 
lution of the 89 GHz data (Figure [JJ. For this analysis, 
a background estimated at 24 counts/pixel was substi- 
tuted for X-ray emission from the pulsar before convolu- 
tion. From the figure, it is clear that the trend of steeper 
emission to the south holds. Although the systematic 
uncertainties are different, this method again implies a 
break in the spectrum of the northern lobe, and allows 
the southern lobe to have a simple power law spectrum 
from radio to X-rays. 

The simplest explanation for the difference in spec- 
trum to the north and south of the pulsar is that we are 
seeing both PWN and shell emission superimposed. In 
this scenario, a majority of the steep-spectrum emission 
seen at 1.4 GHz is attributed to the shell, seen near the 
pulsar due to projection. Other 1.4-GHz emission pre- 
sumably from the shell may be seen at to the east of the 

6 In the 89 GHz image, there are only 3 beam areas within the 
northern region. 
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Fig. 5.— Image s of G29.7-0.3 at (left to right) 1.4 GHz IHelfand et alJI2003blh 89 GHz, and 0.5-10 keV {Chandra archive; as presented 
bv lHelfand et al.ll2003bl) . showing the regions over which emission was integrated in the southern and northern lobes. 



TABLE 2 

Flux densities of the lobes around the pulsar in G29.7— 0.3 



Frequency 
(GHz) 



N. lobe 



Flux Density (uncertainty) 

(mJy) 



y. lobe 



1.4 

89 

1.2 x 10 s (0.5 keV) 



170(35) 
34(5) 
0.00675(5) 



230(35) 
31(5) 
0.00303(3) 



Note. — Datasets for the measurements at 1.4 GHz and 0.5 keV are as referenced in 
Figure [51 
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Fig. 6. — Simple "straight-line" spectra of the northern and 
southern lobes of the PWN in G29.7— 0.3, derived from the data 
in Table 



PWN in Figure E| The extent of the true PWN is then 
seen most easily in X-rays, where the southern portion 
is much fainter. 

The asymmetry in the position of the pulsar with re- 
spect to the northern lobe (the implied true PWN) may 
be due to the passage of a reverse shock from the inter- 
action of the supernova with the surrounding medium. 
An asymmetric medium will cause an asymmetric re- 
verse shock. This phenomenon has been used to account 
for the PWN in the Vela supernova remnant, which is 
largely to one side of the presumed progenitor position 



(Blon din. Chevalier, fc Friersonll2001f) . However, the ex- 
pected 10 4 years until the reverse shock is much larger 
than the age implied by the pulsar's rat e of spin down, 
which falls in the range ~ 700 -1700 years l|Gotthelf et alJ 
l2000HMereghetti et al.ll2002|) . Thus in the absence of an 
extremely dense ISM the passage of a reverse shock does 
not seem a likely explanation for the asymmetry. 

An alternative possibility is that the pulsar may have 
traveled from the apparent center of the PWN since the 
explosion. Assuming that the pulsar was born at the 
brightest part o f the northern lobe, the kinematic dis- 
tance of 21 kpc (Becker & Hclfand 1984) and age range 
above imply a transverse velocity of 400-1000 km s — 1 . 
The upper end of this range is higher than the trans- 
verse velocity of almost all known pulsars, so improved 
pulsar timing (and eventually proper motion measure- 
ments) should allow us to evaluate this scenario in the 
future, although we note that this pulsar has so far been 
detected only in X-rays. 

A final possibility is that the asymmetry in the nebula 
is intrinsic to the pulsar wind. In this case the southern 
lobe could be unassociated shell emission, as suggested 
above, or a part of the PWN with a different spectral 
shape. 

4.4. Origin of the Spectral Breaks 

In both the PWNe in G16. 7+0.1 and G29. 7-0.3, there 
is strong evidence for a low frequency spectral break. 
For G29.7— 0.3, this break may, as mentioned above, oc- 
cur only in the northern lobe of the PWN. Both rem- 
nants are relatively youthful. The conventional spectral 
break due to synchrotron losses in a field B mG occurs 



at age — 405 



-1.5, ,-0.5 



kyr at frequency v GHz. Thus 
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for G16. 7+0.1 we expect a field of ~ 2.5 mG, while for 
G29.7— 0.3 we expect > 1.5 mG. These values are sub- 
stantially higher than inferred in other remnants from 
equipartition arguments, or mo re directly in the Crab 
Nebula (e.g. iHester et al.l ll996). and seem to indicate 
that standard synchrotron losses alone are not the origin 
of the breaks. 

We are left with the conclusion t hat perhaps these ar e 
indeed remnants of a "second kind" fWoltier et al.1 1997'). 
in which intrinsic processes result in a spectrum differ- 
ent from that of the Crab. One possibility is that the 
pulsar 's excitation has decrea sed substantially at some 
point l)Green fc Scheuerll992l). with the break m oving to 
lower frequencies with time ijWoltier et alJll997t) . Alter- 
natively the pulsar injection spectra may not be a simple 
power law. 

5. CONCLUSIONS 

We have carried out observations of three pulsar wind 
nebulae at millimeter wavelengths with the Berkeley- 
Illinois-Maryland Association Array. One of these 
PWNe, in the SNR Gil. 2-0. 3, was not adequately im- 
aged to allow quantitative analysis. We find the break 
frequency of the PWN in G16. 7+0.1 to be at about 26 
GHz, and see an unusual upturn in the spectrum be- 
low or in the X-ray band. In only one other PWN, that 
surrounding PSR B1757— 24, has a similar effect been 
implied. 



The bulk spectrum of the PWN in G29.7-0.3 indicates 
a likely spectral break at 55 GHz, but our observations 
and previous data reveal a marked variation in spectrum 
across the PWN. The spectrum is steeper to the south 
of the pulsar. It may be that the southern emission does 
not represent the continuing effect of the pulsar, but cor- 
responds to a part of the SNR shell seen in projection. 
However, we cannot account for the one-sided PWN un- 
der this scenario. 

The number of PWNe for which an analysis like 
this may be made is limited by the large size of the 
PWNe compared to the spatial scales on which mil- 
limeter interferometers are sensitive. However, the sen- 
sitivity and resolution of next generation instruments, 
such as ALMA, will allow these investigations to be 
conducted on much larger samples in other galaxies. 
Meanwhile, improved single-dish mapping techniques 
and sensitive wide-field imaging interferometers such 
as CARMA could enhance the high-frequency data on 
Galactic PWNe. 
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